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of predictions for two-time response and correlation functions. Applications to non- 
equilibrium critical dynamics of several systems, with a fully disordered initial state and 
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Andersen model and the Ising spin glass are discussed. The critical contact process 
and the parity-conserving non-equilibrium kinetic Ising model are also considered. 
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The analysis of the collective behaviour of many-body systems is greatly helped 
in situations where some scale-invariance allows an efficient description through field- 
theoretical methods. A necessary requirement for the application of these is the 
possibility to identify the physical observables typically defined in terms of a lattice 
model, e.g. a r for the order-parameter at the site r, with a continuum field <p(r) (called 
a scaling operator pQ) with well-defined scaling properties <p(r) = b~ x cf)(r/b). In other 
words, one generally expects that the correspondence (a is the lattice constant) 

a r - a-^(r) (1) 

can be defined in equilibrium systems or more generally steady-states of non-equilibrium 
systems, see e.g. 0IUIS1- In addition, in equilibrium systems one expects the same sort 
of relationship to hold true where <f)(r) is now a primary scaling operator of a conformal 
field-theory and allows space-dependent rescaling factors b = b(r) PQ. 

In this letter, we reconsider this correspondence for systems with dynamical scaling 
and far from equilibrium, as it occurs for example in ageing phenomena. Concrete 
examples are phase-ordering kinetics or non-equilibrium critical dynamics, see [H El IE] 
for reviews. Among the main quantities of interest are the two-time autocorrelation 
function C(t, s) and the autoresponse function R(t, s) 

s- b fc(t/s) 

0(t,r)0( S ,r)) = 8- l -*f R {t/8) (2) 

where <fi is the response field in the Janssen-de Dominicis formalism 0|S], a and b are 
ageing exponents and fc and fn are scaling functions such that fc,n(y) ~ y- x c,n/ z f or 
1/ > 1. These scaling forms are only valid in the scaling regime where t, s — > oo and 
y = t/s > 1 fixed. Until recently, the scaling (J2J) has only been studied for systems with 
a fully disordered initial state with mean initial magnetization mo = (4>(0, r)) = 0. The 
study of the effects of a non- vanishing initial magnetization on the ageing behaviour is 
only beginning [§1 EH E] • We stress that in the kind of system under consideration 
invariance under time-translations is broken. In an attempt to try to derive the form 
of the scaling functions in a model-independent way it has been argued ^2] that the 
scaling operators <fi and <fi should transform covariantly under a larger group than 
mere dynamical scale-transformations. If such an invariance exists, one may call it a 
local scale-invariance (LSI).] The infinitesimal generators of local scale-invariance read 

da mad! 

X = -tdt-- , X 1 = -1?d t --(x + Z)t (3) 

where for simplicity we have suppressed the terms acting on the space coordinates 
which are not important for what follows. We have also not written down the further 

f All existing tests of LSI have been performed for m = 0. 
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generators of LSI which do not modify the time t but only act on the space coordinates 
r, and the absence of any scaling of m means that we are restricting ourselves to 
the case mo = throughout. Here x is the scaling dimension of the scaling operator 
(f>(t, r) = b~ x / z (f)(t/b z , r/b) where z is the dynamical exponent and £ is a constant. It is 
the purpose of this letter to clarify the meaning of this constant £. 

Motivated by the analogy with two-dimensional conformal invariance, we generalize 
the dilatation generator X and the generator X\ of 'special' transformations as follows 
to all n > 

X n = -t n+1 d t - -(n + l)t n - ^-nt n (4) 

such that the commutator [X n , X m ] = (n — m)X n+m holds for all n,m E Nq (with the 
convention G No) 4 Next, the global form of these transformations reads as follows. If 
t = P(t') such that /3(0) = 0, then <f>(t) transforms as 




<j>{t) = (3{t'r X,Z (5) 



where again the space-dependence of <fi was suppressed. The infinitesimal generators 
X n are recovered for (3{t) = t + et n+1 , with |e| 1. From this, it is clear that is not 
transforming as an usual primary scaling operator. But if one defines <&(£) := t~ 2 ^ z (j)(t) 
the scaling operator $(i) becomes a conventional primary scaling operator of LSI, viz. 

$(t) = (3(t')-( x+2 ^' z &(f) (6) 

but with a modified scaling dimension x — > x + 2£. In other words, if time-dependent 
observables of lattice models a r {t) can be related to a primary scaling operator $(t) at 
all, it should be via the relation 

a r (t) -> a~ x (j)(t) = a- x t 2 ^ /z $(t) (7) 

rather than by eq. (pQ). Of course, (J2J) is only possible because of the absence of 
time-translation invariance. We emphasize that the scaling of is unusual in that 
under a dilatation t — > b z t the scaling dimension remains x but for more general scale 
transformations a new effective scaling dimension x + 2£ appears. 

As a simple application, consider the two-time autoresponse function. For quasi- 
primary scaling operators $(£) and $(s) with scaling dimensions x and x, respectively, 
local scale-invariance with mo = predicts ($(t) < l ) (s)) = (t/s)^ x ~ x ^ z (t — s ^-( x + x )/ z ^ U p 
to normalization, as shown in [T2]. In view of (J7J), the physical autoresponse function 
rather reads 



| This is the unique semi-infinite extension of the algebra (Xo,Xi) which does not introduce further 
differential operators into X n and is compatible with eq. Ipjjl. 
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Table 1. Values of the exponents a, a' and Xr/z in several non-disordered and a few 
glassy systems which are at a critical point of their stationary state. If a numerical 
result is quoted without an error bar it is taken form the literature, otherwise the 
numbers in brackets give our estimate of the uncertainty in the last digit (s). nekim 
stands for non-equilibrium kinetic Ising model with conserved parity and FA stands 
for the Frederikson-Andersen model. The methods of calculation of the two-time 
autoresponse are D: direct space, P: momentum space, A: alternating external field; E 
refers to an exact solution and N to a numerical study. 



model 


a 


a' — a 




Method 


Ref. 


OJK-model 


(d-l)/2 


-1/2 


d/A 


D,E 


17, 18, 14 


ID Ising 





-1/2 


1/2 


D,E 


nap 


2D Ising 


0.115 


-0.187(20) 


0.732(5) 


P,N 


22 


3D Ising 


0.506 


-0.022(5) 


1.36 


P,N 


22 


ID contact process 


-0.681 


+0.270(10) 


1.76(5) 


D,N 


26, 27J 


ID NEKIM 


-0.430 


-0.09 


0.56 


D,N 


[28 


FA, d > 2 


l + d/2 


-2 


2 + d/2 


P,E 




FA, d = 1 


1 


-3/2 


2 


P,E 


[2011211 


3D Ising spin glass 


0.060(4) 


-0.76(3) 


0.38(2) 


A,N [14 



R(t,s) = (<P(t)<P(s)) = (t 2 ^$(t)^>$( S )} 



[2£+x-x)/z \ ~{x+x+2(,+2£)/z 



| - I I- - ! 

l+a'—Xn/z / ^ \ — 1— a' 



= ^°(;J [s- 1 ) < 8 > 

(up to normalization) and the effective scaling dimensions of $(£) and $(s) are read 
off from eq. (@) to be now x + 2£ and x + 2£, respectively. In the last line, we have 
reintroduced the standard exponents a, a' and Xr and hence reproduce the result quoted 
in [T^j. Early discussions of local scale-invariance had assumed a' = a from the outset. 
In the appendix, we discuss the scaling form of the autocorrelator C(t, s) in those cases 
where z — 2. 

It appears that the more general correspondence (J7J) and consequently the response 
© with a' ^ a actually occurs in non-equilibrium critical dynamics, as we shall now 
illustrate in a few examples. We stress that in the models considered here (with the 
only exception of the contact process) we always use a fully disordered initial state 
with a vanishing initial magnetization mo = 0.§ In table ^ we collect results on the 
exponents a, a' and Xr/z in some models with a critical stationary state and where 

§ From LSI, it is then easy to see that the time-dependent magnetization m(t) = mo = for all times, 
in agreement with the Monte Carlo and the exact results. On the other hand, if initially tuq > 0, one 
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Figure 1. Intermediate susceptibility xim(t, s) in momentum space in the (a) 2D and 
(b) 31? critical Ising model, for several values of the waiting time s. The full curve is 
the LSI prediction eq. Hl()lll|l with the exponents as listed in tabled The dashed line 
corresponds to the case a' = a. 



a' a. || In several cases, these exponents can be read off from the exact solution, 
i.e., for the magnetic response in the OJK-model [T7| H8] and the ID Glauber-Ising 
model at zero temperature ^H] or else the energy response in the zero-temperature 
Frederikson- Andersen model 



Another interesting test case is provided by the critical Ising model in 2D and 3D. 
Indeed, it was pointed out some time ago that the numerical calculation of the two-time 
response function R q (t, s) = / M ddr R(t, s; r)e~ iq r in momentum space provides a more 
sensitive test on the form of its scaling function than in direct space [22] • The methods 
of LSI can be readily adapted to momentum space and the analogue of (JBJ) is, again up 
to normalization and for mo = 

(,\l+a'—X R /z /, \—X—a'+d/z 
s) (j " 1 ) (9) 

Since measurements of autoresponse functions are much affected by statistical noise, 
one often rather studies integrated response functions. Here we consider 

s 



XM(t,s):= duR (t,u)=xos~ a+d/z f x (t/s) (10) 



-a+d/z 

A.** j-i-uv") "7 All" j 

'a/2 

has the regime of short-time dynamics with m(t) ~ t s ^H] before the long-time decay m(t) ~ t~i i l { ~ lJZ ) 
The scaling of two-time observables has been recently discussed in [^l llOlITT] and it was shown that 
for tuq the universal scaling behaviour is different from the one found for mo = 0. An extension of 
LSI to non-equilibrium critical dynamics with non-vanishing initial magnetizations is an open problem 
to which to hope to return elsewhere. 

|| In tabled D > E means that the exact response agrees with (jHJl with the given values of the exponents, 
while p,e means that there is exact agreement with (jHJl. 
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Figure 2. Autoresponse function for the critical ID contact process for several waiting 
times s. The data labelled TM come from the transfer matrix renormalization group 
and MC denotes Hinrichsen's Monte Carlo data [22| - The dashed line corresponds 
to the case a' = a and the full curve gives the LSI prediction eq. (jHJ with the exponents 
as listed in tabled 



which is free from effects which mask the true scaling behaviour in several other variants 
of integrated responses [22J. The scaling function f x (y) follows from LSI, eq. (jHJ): 



and where 2F1 is Gauss' hypergeometric function. In figure ^ we compare simulational 
data [22] with this prediction for both the 2D and 3D critical Ising model with non- 
conserved heat-bath dynamics. It had already been observed before [22] that local scale- 
invariance with the additional assumption a' = a does not agree with the numerical data 
in 2D and only marginally so in 3D and we confirm this finding. However, we also see 
that the data can be perfectly matched by LSI, within the numerical precision, if a and 
a' are allowed to be different. We did check that the integrated TRM response functions 
in direct space as studied in [23] do not change appreciably with a' — a. 

A similar conclusion can also be drawn for the ID critical contact process. It has 
been shown recently that the phenomenology of ageing can also be found in critical 
stochastic processes although these do not satisfy detailed balance and have a non- 
equilibrium steady-state [2HI23I2S]- m figure|2]we compare the numerical data obtained 
directly for R(t, s) either from the LCTMRG [26J or Monte Carlo simulations [2Z|. It is 
satisfying that the data from both methods are consistent with each other in the scaling 
regime, where s and t — s are both large enough. Again, we observe an almost perfect 




(11) 
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agreement with eq. (|SJ), provided a' ^ 

One the other hand, when one looks closer at the region where t/s J$ 1.1, one does 
observe deviations of the data from (jHJ) . In trying to analyze this, recall that non- 
equilibrium critical dynamics is special in the sense that both the ageing regime (where 
t — s ~ O(s)) and the quasistationary regime (where t — s <C s) display dynamical 
scaling with the same length scale L{t) ~ t 1 ^, where z is the equilibrium dynamical 
critical exponent. Hence one usually expects some 'crossover' to occur. In terms of the 
response function, this might be formalized by writing R = IZ^s/t*, (t — s)/r*, s) where 
r* is some reference time scale such that, with (t — s)/r* = 0(1) 

11m d = I ^ q (t - s) ; for s/n -> oo 

\ s-i-Mt/s) ; for s/n = 0(1) 

Since in lattice calculations s is always finite, the 'crossover' can be illustrated by 
studying Q := R(t , s) / R eq (t — s) ~ R(t,s)(t — s) 1+a . As long as LSI still describes 
the data, one expects Q ~ (y — l) a_a ' for y = t/s <J 1 and deviations from it should 
signal the presumed crossover to the quasistationary regime y — > 1 where Q(y) should 
become constant. For the ID critical contact process we find that for y — t/s ^ 1.1, Q(y) 
still obeys scaling for s large enough and that Q changes from Q w 0.3 at y — 1 ~ 0.1 
to Q ~ 0.8 at y — 1 ~ 10 -4 . <3(y) appears to become flatter as y — > 1, but the 
change to a quasistationary behaviour could not yet be observed, in spite of large 
waiting times s > 860000, before strong finite-time effects set in at t — s = 0(1). 
In comparison, unpublished data for the 2D Ising model ^Zj show convergence towards 
Qiy) ~ (y — l) 0187 as s increases before finite-time effects destroy scaling. We conclude 
that LSI does accurately describe the data as long as t/s is large enough such that the 
effects of the 'crossover' are not yet noticeable. A quantitative analysis of data from the 
region t/s J$ 1.1 would require a precise theory of the 'crossover' between the ageing 
regime and the region t-s <C s, the role of finite-time effects and the influence of initial 
conditions, e.g. different initial fillings of the lattice. 

Very recently, a similar test was carried out in a ID kinetic Ising model with 
competing Glauber and Kawasaki dynamics [2B1- The stationary state is therefore not 
an equilibrium state. This was the first time that LSI was tested and confirmed for a 
dynamics where the parity of the total spin is conserved by the dynamics. 

Finally, we recall that studying the scaling behaviour of an alternating susceptibility 
gives yet another direct access to the exponent a' — a. This was applied to the critical 
3D Ising spin glass [T3j, with a binary distribution of the couplings = ±J. 

In summary, we have reconsidered the way how observables defined in non- 
equilibrium lattice models might be related to (quasi-) primary scaling operators of 
field-theory. Our result eq. © points to a so far overlooked subtlety which might be of 

% Hinrichsen quotes Xr/z « 1.75 and 1 + a' « 0.59 [22] in good agreement with our estimates. The 
contact process is the only known example where a' — a > 0. 
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relevance in the discussion of the functional form of non-equilibrium scaling functions, 
for example in ageing phenomena. It remains to be seen how general the phenomenon 
for which we have presented evidence really is. + The results on R(t, s) as collected 
in table ^ for the non-equilibrium dynamics of some models with m = appear to 
be compatible with the predictions eqs. ()8J9j) of local scale-invariance, provided cross- 
over effects to non-ageing regimes are negligible. The multitude of examples in tabled 
suggests that rather being a kind of exotic exception (a belief implicit in O EH EI), 
the case a' ^ a might turn out to be the generic situation. Having seen that the 
same mechanism also explains the exact autocorrelator of the ID Glauber-Ising model 
indicates that the correspondence (J2|) should be more than just a patching-up of data 
for the autoresponse function. 

What does this mean for the existence of local scale-invariance in non-equilibrium 
dynamics, with tuq = ? In a few exactly solved systems (where the dynamical exponent 
z = 2) we have found exact agreement and in several models as generic as kinetic Ising 
models or the contact process eqs. (|8I9|) describe the data very well for t/s not too 
small. It is remarkable that the two-time autocorrelations and autoresponses of models 
as physically different as those included in table ^ (and several further ones with a = a' 
which we did not include) can be described in terms of a single theoretical idea. On the 
other hand, field-theoretical studies of the critical O(n) model in both 4 — e dimensions 
[5[l29] and in 2 + e dimensions [Tl], although they agree with LSI at the lowest orders in 
e, continue to find discrepancies with either (JHJ) or © at some higher order. However, 
non-equilibrium field-theory presently only yields explicit results for the first few terms 
of the e-expansion series. When one truncates this series to an e-dependent sum, the 
resulting numerical values for the scaling functions are still far from the numerical data.* 
But since we have shown that LSI reproduces the known exact results of both R(t, s) and 
C(t, s) of the ID Ising model it might be too simplistic to argue that LSI could at best 
describe gaussian fluctuations. A better understanding of the dynamical symmetries of 
non-equilibrium critical dynamics remains a challenging problem. 

Appendix. Two-time autocorrelations for z = 2 

If the dynamical exponent z = 2, local scale-invariance reduces to Schrodinger- 
invariance. We have already described in the past [13] how two-time autocorrelation 
functions can be calculated in the case £ = and we now wish to extend that treatment 
to the more general correspondence (J2J). We consider a Langevin equation of the form 
d t (f) = —Djl — Dv(t)<f> + 1] where 7i is the hamiltonian, D the diffusion constant, the 
gaussian noise r\ has zero mean and variance {rj(t,r)r](s,r')) = 2DT5(t — s)8(r — r') 

+ It is not inconceivable that analogues might exist in equilibrium critical phenomena, for instance 
when spatial translation-invariance is broken by disorder or boundaries. 

* The second-order calculation in 4 — e dimensions for n — 1 is a little closer to the numerical data 
than LSI with a' = a [22|. 
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and T is the bath temperature. The potential v(t) acts as a Langrange multiplier which 
can be used to describe explicitly the breaking of time-translational invariance. Here 
we restrict to situations where 



k(t) := exp 



D / duv(u) 



o 



(Al) 



Then is has been shown jTHI that for systems at criticality 

C{t, s) = ) = DT cj dR (<f>(t, y)(j){s, y)4>\u, r + y] 



DT r dudR' 



k{u) 



-K ( Q>(t,s,u;R) 



(A2) 



/0\ 

where 1Z Q is the well-known three-point response function for v(t) — which can be 
found from its Schrodinger-covariance and reads [3*U] 



(t, s, u; r) = IZ^' (t, s, u) exp 



(3), 



M t + s - 2u 



(t — u)(s — u) 



2 (s -u)(t-u) 

nf\t, s, u) = Q{t - u)Q{s -u){t- u)~* 2 {s - u)-^{t - s)~ x+ ^ 

where \I/ is an undetermined scaling function and the causality conditions t > u, s > u 
are noted explicitly In writing this, we have dropped a term coming from the 
correlations in the initial state which merely produces finite-time corrections to the 
leading scaling behaviour, see jH EJ H3] ■ 

We now generalize this to the primary scaling operators according to ((7j). The 
operator $ has the scaling dimension x + 2£ and the composite scaling operator $ 2 
has the scaling dimension 2x2 + We then obtain for the physical autocorrelation 
function, up to normalization and with t > s 

C(t,s) = (ts)t [dudR ($(t,y)$(s,y)$ 2 (u,R + y)) vft 2 



du 



dR exp 

^l+d/2-x-x 2 



Mt + s - 2m o 
— R 



^^-u 2 ^ [(t -u)(s- 
* (R 2 ) 



s J 



dvv^ 2 ' F 



— V 



x 2 +2&-x-2Z-d/2 



d/2-x 2 -2& 



t/s+ l-2v 
t/s - 1 



(A3) 



and where the function ^ is defined by the integral over R. By comparison with the 
standard scaling from for C(t, s), we read off b = x+x 2 — 1— d/2 and Ac = 2{x— / r )+2^.f f 

ft For bosonic free fields, one would have X2 = x and £2 = £■ 

ff A similar calculation for the autoresponse function gives, up to normalization, 



R(t, s) = s-( x+ *V 2 (t/s)S +r (t/s - l)- x ~ 2( 5, 
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Furthermore, since 1 + a' = x + 2£, it turns out that the form of the scaling function 
fc{y) is described by just one more parameter fi := £ + £ 2 and we finally have 

/,\ l+a'-A c /2 /, \ b-2a'-l+2fJ, 



X 



^ v v X C +2fl-2-2a' 







--v) 1- 



and we have also reintroduced a normalization constant Cq. This should hold for simple 
(non-glassy) magnets with z = 2 and in situations where the initial correlations have 
no effect on the leading scaling behaviour; of course the scaling limit s —>■ oo and 
t/s = y > 1 fixed is understood. 

As an illustration, we consider the ID Glauber-Ising model. At T = 0, the exact 
two-time autocorrelation function is ^H] 



C(t,s) = — arctan I I (A5) 



7T 



t/s-l 



This holds true not only for the usually considered short-ranged initial conditions but 
also for long-ranged initial spin-spin correlations (ov(0)<To(0)) ~ r~ u with v > (for 
v = an analogous result holds for the connected autocorrelator) [T^]. In addition, the 
exponents a, a' and Xr are independent of v. 

In previous work ^3] , we have already explained the form of the exact autoresponse 
function R(t, s) in terms of the correspondence eq. (J7J) (see table HJ) but we had to leave 
open the analogous question for C(t, s). In order to account for (|A5J) . we remark that 
for t = s, the autocorrelator should not be singular. This requires 

V( W ) = w b-2a>-l+2» for w > x ^ 

The most simple way to realize this is to require that ()A6|1 holds for all values of w. 
This kind of assumption was already seen to become exact in models described by an 
underlying bosonic free field-theory ^Hj • Recalling from table Q that b = a = and 
Xc = 1 and assuming (jA6|) to hold for all w, we obtain 



C(t,s) « C [ d^ 2 " 
Jo 



-2^-1/2 / ^ \ 2/j 



1 - 2v) (A7) 



Because the exact result (jA5|) is independent of the initial correlations, the comparison 
with the expression (jA7|) derived from the thermal noise is justified. The exact Glauber- 
Ising result (jA5|) is recovered from (|A7|) for /i = —1/4 and Cq = \/2/%. 

which reproduces again eq. JHJ), hence A^ = 2(x — F) + 2£ = Ac as expected 0] for non- 
disordered systems without long-range initial correlations. In particular, for critical systems with 
a = b the equality Ac = A# implies that there is a finite limit fluctuation-dissipation ratio = 
Um (t/s) ^ 00 R(t, s)/{T c dC{t, 8)188), see 0. 
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This is the first example of an exactly solved model with a' ^ a where the scaling 
of both the autoresponse and of the autocorrelation functions can be explained in terms 
of LSI. 
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